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Abstract — Underwater  acoustic  (UWA)  channels  are  wideband 
in  nature  due  to  the  small  ratio  of  the  carrier  frequency  to  the 
signal  bandwidth,  which  introduces  frequency-dependent  Doppler 
shifts.  In  this  paper,  we  treat  the  channel  as  having  a  common 
Doppler  scaling  factor  on  all  propagation  paths,  and  propose  a 
two-step  approach  to  mitigating  the  Doppler  effect:  1)  nonuniform 
Doppler  compensation  via  resampling  that  converts  a  “wideband” 
problem  into  a  “narrowband”  problem  and  2)  high-resolution 
uniform  compensation  of  the  residual  Doppler.  We  focus  on 
zero-padded  orthogonal  frequency-division  multiplexing  (OFDM) 
to  minimize  the  transmission  power.  Null  subcarriers  are  used  to 
facilitate  Doppler  compensation,  and  pilot  subcarriers  are  used 
for  channel  estimation.  The  receiver  is  based  on  block-by-block 
processing,  and  does  not  rely  on  channel  dependence  across 
OFDM  blocks;  thus,  it  is  suitable  for  fast- varying  UWA  channels. 
The  data  from  two  shallow-water  experiments  near  Woods  Hole, 
MA,  are  used  to  demonstrate  the  receiver  performance.  Excellent 
performance  results  are  obtained  even  when  the  transmitter  and 
the  receiver  are  moving  at  a  relative  speed  of  up  to  10  kn,  at  which 
the  Doppler  shifts  are  greater  than  the  OFDM  subcarrier  spacing. 
These  results  suggest  that  OFDM  is  a  viable  option  for  high-rate 
communications  over  wideband  UWA  channels  with  nonuniform 
Doppler  shifts. 
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1.  Introduction 

ULTICARRIER  modulation  in  the  form  of  orthogonal 
frequency-division  multiplexing  (OFDM)  has  prevailed 
in  recent  broadband  wireless  radio  applications  due  to  the  low 
complexity  of  receivers  required  to  deal  with  highly  disper¬ 
sive  channels  [2],  [3].  This  fact  motivates  the  use  of  OFDM 
in  underwater  environments.  Earlier  works  on  OFDM  focus 
mostly  on  conceptual  system  analysis  and  simulation-based 
studies  [4]-[7],  while  experimental  results  are  extremely  scarce 
[8]-[12].  Recent  investigations  on  underwater  OFDM  commu¬ 
nication  include  [13]  on  noncoherent  OFDM  based  on  ON-OFF 
keying,  [14]  on  a  low-complexity  adaptive  OFDM  receiver,  and 
[15]  on  a  pilot- tone-based  block-by-block  receiver. 

In  this  paper,  we  investigate  the  use  of  zero-padded  OFDM 
(ZP-OFDM)  [2],  [16]  for  UWA  communications.  Zero  padding 
is  used  instead  of  cyclic  prefix  to  save  the  transmission  power 
spent  on  the  guard  interval.  The  performance  of  a  conventional 
ZP-OFDM  receiver  is  severely  limited  by  the  intercarrier  in¬ 
terference  (ICI)  induced  by  fast  channel  variations  within  each 
OFDM  symbol.  Furthermore,  the  UWA  channel  is  wideband 
in  nature  due  to  the  small  ratio  of  the  carrier  frequency  to  the 
signal  bandwidth.  The  resulting  frequency-dependent  Doppler 
shifts  render  existing  ICI  reduction  techniques  ineffective.  We 
treat  the  channel  as  having  a  common  Doppler  scaling  factor  on 
all  propagation  paths,  and  propose  a  two-step  approach  to  mit¬ 
igating  the  frequency-dependent  Doppler  shifts:  1)  nonuniform 
Doppler  compensation  via  resampling,  which  converts  a  “wide¬ 
band”  problem  into  a  “narrowband”  one  and  2)  high-resolution 
uniform  compensation  of  the  residual  Doppler  for  best  ICI  re¬ 
duction. 

The  proposed  practical  receiver  algorithms  rely  on  the  pre¬ 
amble  and  postamble  of  a  packet  consisting  of  multiple  OFDM 
blocks  to  estimate  the  resampling  factor,  the  null  subcarriers 
to  facilitate  high-resolution  residual  Doppler  compensation, 
and  the  pilot  subcarriers  for  channel  estimation.  The  receiver 
is  based  on  block-by-block  processing,  and  does  not  rely  on 
channel  coherence  across  OFDM  blocks;  thus,  it  is  suitable  for 
fast- varying  underwater  acoustic  (UWA)  channels.  To  verify 
our  approach,  two  experiments  were  conducted  in  shallow 
water:  one  in  the  Woods  Hole  Harbor,  MA,  on  December  1, 
2006,  and  the  other  in  Buzzards  Bay,  MA,  on  December  15, 
2006.  Over  a  bandwidth  of  12  kHz,  the  data  rates  are  7.0,  8.6, 
9.7  kb/s  with  quaternary  phase- shift  keying  (QPSK)  modula¬ 
tion  and  rate  2/3  convolutional  coding,  when  the  numbers  of 
subcarriers  are  512,  1024,  and  2048,  respectively.  Excellent 
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performance  is  achieved  for  the  latter  experiment,  while  rea¬ 
sonable  performance  is  achieved  for  the  former  experiment 
whose  channel  has  a  delay  spread  much  longer  than  the  guard 
interval.  The  receiver  performs  successfully  even  at  a  relative 
speed  of  up  to  10  kn,  resulting  in  Doppler  shifts  that  are  greater 
than  the  OFDM  subcarrier  spacing.  These  results  suggest  that 
OFDM  is  a  viable  option  for  high-rate  UWA  communications 
over  UWA  channels. 

The  rest  of  this  paper  is  organized  as  follows.  In  Section  II, 
the  performance  of  a  conventional  OFDM  receiver  is  analyzed. 
In  Section  III,  a  two-step  approach  to  mitigating  the  Doppler 
shifts  is  proposed,  and  the  practical  receiver  algorithms  are  spec¬ 
ified.  In  Sections  IV  and  V,  the  receiver  performance  is  reported. 
Section  VI  contains  the  conclusions. 

II.  Zero-Padded  OFDM  eor  UWA  Channels 

Let  T  denote  the  OFDM  symbol  duration  and  Tg  the  guard 
interval.  The  total  OFDM  block  duration  is  T'  =  T  -h  T^.  The 
frequency  spacing  is  A/  =  1/T.  The  kih  subcarrier  is  at  the 
frequency 


is  the  case  in  our  experiments,  this  assumption  seems  to 
be  justified. 

A2)  The  path  delays  r^,  the  gains  and  the  Doppler 
scaling  factor  a  are  constant  over  the  block  duration  T' . 
The  OFDM  block  durations  are  T  =  42.67,  85.33,  and 
170.67  ms  in  our  experiments  when  the  numbers  of  sub¬ 
carriers  are  512,  1024,  and  2048,  respectively.  Assump¬ 
tion  A2)  is  reasonable  within  these  durations,  as  the 
channel  coherence  time  is  usually  on  the  order  of  sec¬ 
onds. 

The  received  signal  in  passband  is  then 

y{t)  =  Re  I Z] 

I  p  IkeSA 

X  +n(i)  (5) 

where  h{t)  is  the  additive  noise.  The  baseband  version  y(t)  of 
the  received  signal  satisfies  y{t)  =  and  can  be 

written  as 


fk  =  fc  +  kAf,  k  =  -K/2, . . .  ,iT/2  -  1  (1) 


where  fc  is  the  carrier  frequency  and  K  subcarriers  are  used  so 
that  the  bandwidth  is  B  =  KAf. 

Let  us  consider  one  ZP-OFDM  block.  Let  d[k]  denote  the  in¬ 
formation  symbol  to  be  transmitted  on  the  k\h  subcarrier.  The 
nonoverlapping  sets  of  active  subcarriers  Sa  and  null  subcar¬ 
riers  Sn  satisfy  Sa  U  Sn  =  {— iT/2, . . . ,  iT/2  —  1}.  The  trans¬ 
mitted  signal  in  passband  is  then  given  by 


s{t)  =  Re 


.keSA 


fct 


t  G  [0,T  +  T^: 


(2) 


where  g{t)  describes  the  zero-padding  operation,  i.e.,  g{t)  = 
^  [^^T]  and  g{t)  =  0,  otherwise. 

We  consider  a  multipath  underwater  channel  that  has  the  im¬ 
pulse  response 


c{T,t)  =  'Y^Ap{t)8{T-Tp{t))  (3) 

p 


yit)  = 

k^SA 


^j27vkAft 


^j27Tafkt 


X 


-  P 


+  n{t)  (6) 


where  n{t)  is  the  additive  noise  in  hasehand.  Based  on  the  ex- 
pression  in  (6),  we  observe  the  following  two  effects. 

1)  The  signal  from  each  path  is  scaled  in  duration,  from  T 
to  T/ ( 1  -h  a). 

2)  Each  subcarrier  experiences  a  Doppler-induced  fre¬ 
quency  shift  which  depends  on  the  frequency 

of  the  subcarrier.  Since  the  bandwidth  of  the  OFDM 
signal  is  comparable  to  the  center  frequency,  the 
Doppler-induced  frequency  shifts  on  different  OFDM 
subcarriers  differ  considerably;  i.e.,  the  narrowband 
assumption  does  not  hold. 

The  frequency-dependent  Doppler  shifts  introduce  strong 
intercarrier  interference  if  an  effective  Doppler  compensation 
scheme  is  not  performed  before  the  OFDM  demodulation. 


where  Ap{t)  is  the  path  amplitude  and  rp{t)  is  the  time- varying 
path  delay.  To  develop  our  receiver  algorithms,  we  adopt  the 
following  assumptions. 

Al)  All  paths  have  a  similar  Doppler  scaling  factor  a  such 
that 

Tp{t)  iHTp-  at.  (4) 

In  general,  different  paths  could  have  different  Doppler 
scaling  factors.  The  method  proposed  in  this  paper  is 
based  on  the  assumption  that  all  the  paths  have  the  same 
Doppler  scaling  factor.  When  this  is  not  the  case,  part  of 
useful  signals  are  treated  as  additive  noise,  which  could 
increase  the  overall  noise  variance  considerably.  How¬ 
ever,  we  find  that  as  long  as  the  dominant  Doppler  shift 
is  caused  by  the  direct  transmitter/receiver  motion,  as  it 


III.  Receiver  Design 

First,  we  present  in  Section  II- A  the  technical  approach  to 
mitigating  the  frequency-dependent  Doppler  shifts,  and  then 
specify  in  Section  III-B  practical  receiver  algorithms  that  we 
apply  to  the  experimental  data. 

A.  Two-Step  Approach  to  Mitigating  the  Doppler  Effect 

We  propose  a  two-step  approach  to  mitigating  the  frequency- 
dependent  Doppler  shifts  due  to  fast- varying  UWA  channels  as 
follows. 

1)  Nonuniform  Doppler  compensation  via  resampling.  This 
step  converts  a  “wideband”  problem  into  a  “narrowband” 
problem. 

2)  High-resolution  uniform  compensation  of  residual 
Doppler.  This  step  fine  tunes  the  residual  Doppler  shift 
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Fig.  I.  Detailed  receiver  diagram  on  one  receive  element. 


corresponding  to  the  “narrowband”  model  for  best  ICI 
reduction. 

The  resampling  methodology  has  been  shown  effective  to 
handle  the  time-scale  change  in  underwater  communications; 
see,  e.g.,  [17]  and  [18].  Resampling  can  be  performed  either  in 
passband  or  in  baseband.  For  convenience,  let  us  present  these 
steps  using  passband  signals.  In  the  first  step,  we  resample  the 
received  waveform  y{t)  using  a  resampling  factor  h 

*W  =  1'(tTs)^ 

Resampling  has  two  effects:  1)  it  rescales  the  waveform  and  2)  it 
introduces  a  frequency-dependent  Doppler  compensation.  With 
y{t)  from  (5)  and  z{t)  =  the  baseband  signal 

z{t)  is 


Compensating  for  the  CFO  in  z{t),  we  obtain 

keSA  I  V  p 


(11) 


where  the  subcarriers  stay  orthogonal.  On  the  output  of  the  de¬ 
modulator  in  the  mth  subchannel,  we  have  [2],  [16] 

Tg+T 

Zra  =  ^  j 

0 

~  C{fm)d[m]  +  Vm  (12) 


z{t)  = 


X  <  d\k\e^ 
keSA  I 


,j27TkAf^t 


X 


E4 


-j27rfkTp 


L  P 


1  a 


■v{t) 


(8) 


where  v(t)  is  the  additive  noise.  The  target  is  to  make  (l-fa)/ 
(1  -h  6)  as  close  to  one  as  possible.  With  this  in  mind,  we  have 


z{t)  ^ 


X 


E 


J27rkAft 


k^SA 

+  v{t). 


(9) 


The  residual  Doppler  effect  can  be  viewed  as  the  same  for  all 
subcarriers.  Hence,  a  wideband  OFDM  system  is  converted  into 
a  narrowband  OFDM  system  with  a  frequency-independent 
Doppler  shift 


e  = 


a  —  b 

m 


/c 


(10) 


In  radio  applications,  a  carrier  frequency  offset  (CFO)  between 
the  transmitter  and  the  receiver  leads  to  an  expression  of  the 
received  signal  in  the  form  (9)  [19],  [20].  For  this  reason,  we 
call  the  term  e  in  (10)  as  CFO  when  a  narrowband  model  is 
concerned. 


where  C{f)  :=  and  Vm  is  the  resulting  noise. 

Hence,  ICI-free  reception  is  approximately  achieved.  Rescaling 
and  phase  rotation  of  the  received  signal  thus  restores  the  or¬ 
thogonality  of  the  subcarriers  of  ZP-OFDM.  The  correlation  in 
(12)  can  be  performed  by  overlap-adding  of  the  received  signal, 
followed  by  fast  Fourier  transform  (FFT)  processing  [2],  [16]. 

In  practice,  the  scale  factor  h  and  the  CFO  e  need  to  be  de¬ 
termined  from  the  received  data.  They  can  be  estimated  either 
separately  or  jointly.  Note  that  each  estimate  of  h  will  be  associ¬ 
ated  with  a  resampling  operation,  which  is  costly.  It  is  desirable 
to  limit  the  number  of  resampling  operations  to  as  few  as  pos¬ 
sible.  At  the  same  time,  high-resolution  algorithms  are  needed 
to  fine  tune  the  CFO  term  e  for  the  best  ICI  reduction. 

Next,  we  specify  the  practical  algorithms  that  we  apply  to  the 
experimental  data. 

B.  Practical  Receiver  Algorithms 

The  received  signal  is  directly  sampled  and  all  processing  is 
performed  on  discrete-time  entries.  Fig.  1  depicts  the  receiver 
processing  for  each  element,  where  BPF,  LPF,  and  VA  stand 
for  bandpass  filtering,  lowpass  filtering,  and  Viterbi  algorithm, 
respectively.  Next,  we  discuss  several  key  steps. 

1 )  Doppler  Scaling  Factor  Estimation:  Coarse  estimation  of 
the  Doppler  scaling  factor  is  based  on  the  preamble  and  the  post¬ 
amble  of  a  data  packet.  (This  idea  was  used  in,  e.g.,  [17]  for 
single-carrier  transmissions.)  The  packet  structure,  containing 
Nh  OFDM  blocks,  is  shown  in  Fig.  2.  By  cross  correlating  the 
received  signal  with  the  known  preamble  and  postamble,  the  re¬ 
ceiver  estimates  the  time  duration  of  a  packet  The  time  du¬ 
ration  of  this  packet  at  the  transmitter  side  is  Ttr-  By  comparing 
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Fig.  2.  Packet  structure. 


Trx  with  Ttx,  the  receiver  infers  how  the  received  signal  has  been 
compressed  or  dilated  by  the  channel 


Ttx 

\  Qj 


(13) 


The  receiver  then  resamples  the  packet  with  a  resampling  factor 
b  =  a  used  in  (7).  We  use  the  polyphase-interpolation-based 
resampling  method  available  in  Matlab. 

2)  CFO  Estimation:  A  CFO  estimate  is  generated  for 
each  OFDM  block  within  a  packet.  We  use  null  subcarriers 
to  facilitate  estimation  of  the  CFO.  We  collect  K  L  sam¬ 
ples  after  resampling  for  each  OFDM  block  into  a  vectori 
z  =  [2^(0), . . . ,  z{K  L  —  1)]^,  assuming  that  the  channel  has 
1/  -h  1  taps  in  discrete  time.  The  channel  length  can  be  inferred 
based  on  the  synchronization  output  of  the  preamble,  and 
its  estimation  does  not  need  to  be  very  accurate.  We  define  a 
(K+L)  X 1  vector  =  [1, 

and  2i  (K  L)  X  (K  L)  diagonal  matrix  r(e)  = 

where  =  T/K 
is  the  time  interval  for  each  sample.  The  energy  of  the  null 
subcarriers  is  used  as  the  cost  function 


J(e)=  y]  |Cr«(e)z 


(14) 


If  the  receiver  compensates  the  data  samples  with  the  correct 
CFO,  the  null  subcarriers  will  not  see  the  ICI  spilled  over  from 
neighboring  data  subcarriers.  Hence,  an  estimate  of  e  can  be 
found  through 


e  =  argminJ(e)  (15) 

which  can  be  solved  via  1-D  search  for  e.  This  high-resolution 
algorithm  corresponds  to  the  MUSIC-like  algorithm  proposed 
in  [19]  for  cyclic-prefixed  OFDM. 

Instead  of  the  1-D  search,  one  can  also  use  the  standard  gra¬ 
dient  method  as  in  [20]  or  a  bisectional  search.  A  coarse-grid 
search  is  needed  to  avoid  local  minima  before  the  gradient 
method  or  the  bisectional  search  is  applied  [21]. 

Remark  1:  The  null  subcarriers  can  also  facilitate  joint  re¬ 
sampling  and  CFO  estimation.  This  approach  corresponds  to  a 
2-D  search:  when  the  scaling  factor  h  and  the  CFO  e  are  cor¬ 
rect,  the  least  signal  spillover  into  null  subcarriers  is  observed. 
However,  the  computational  complexity  is  high  for  a  2-D  search. 
This  algorithm  can  be  used  if  no  coarse  estimate  of  the  Doppler 
scaling  factor  (e.g.,  from  the  pre-  and  postamble  of  a  packet)  is 
available. 

^Bold  upper  case  and  lower  case  letters  denote  matrices  and  column  vectors, 
respectively;  (•)*,  and  (  denote  transpose,  conjugate,  and  Hermitian 
transpose,  respectively. 


3)  Pilot-Tone-Based  Channel  Estimation:  After  resampling 
and  CFO  compensation,  the  ICI  induced  by  CFO  is  greatly  re¬ 
duced.  Due  to  assumption  A2),  we  will  not  consider  the  ICI 
because  of  channel  variations  within  each  OFDM  block.  Note 
that  ICI  analysis  and  suppression  in  the  presence  of  fast- varying 
channels  have  been  treated  extensively  in  the  literature;  see,  e.g., 
the  references  listed  in  [22,  Ch.  19].  Ignoring  ICI,  the  signal  in 
the  mth  subchannel  can  be  represented  as  [c.f.,  (12)] 

Zm  =  f,^r^(e)z  =  H{m)d[m]  +  Vm  (16) 


where  H{m)  =  C{fm)  is  the  channel  frequency  response  at 
the  mth  subcarrier  and  Vm  is  the  additive  noise.  On  a  multipath 
channel,  the  coefficient  H{m)  can  be  related  to  the  equivalent 
discrete-time  baseband  channel  parameterized  by  T  -h  1  com¬ 
plex-valued  coefficients  through 

L 

H{m)  =  (17) 

1=0 


To  estimate  the  channel  frequency  response,  we  use  pilot 
tones  at  subcarrier  indices  pi,...,PiTp;  i-e.,  {d[pi]}iJi  are 
known  to  the  receiver. 

As  long  as  iCp  >  T  -h  1,  we  can  find  the  channel  taps  based 
on  a  least  squares  (LS)  formulation 


"  ' 

'  Vpi  ' 

'd\pi] 

= 

+ 

-^PKp  _ 

ypKp  _ 

^  [PK,]  _ 

•  — Zp  ; — D 


'1  ••• 

~  ho  ~ 

\  g— . 

Q—j^PKpL 

-hL- 

:=V 

:=h 

(18) 


To  minimize  the  complexity,  we  will  adhere  to  the  following 
two  design  rules: 

dl)  the  Kp  pilot  symbols  are  equally  spaced  within  K  sub¬ 
carriers; 

d2)  the  pilot  symbols  are  PSK  signals  with  unit  amplitude. 
Since  the  pilots  are  equispaced,  we  have  that  V  =  KplL^i 
[23],  and  since  they  are  of  unit  amplitude,  we  have  that 
=  Ik,  .  Therefore,  the  LS  solution  for  (18)  simplifies  to 

hLS  =  Ey^B^Zp.  (19) 

IVp 

This  solution  does  not  involve  matrix  inversion,  and  can  be  im¬ 
plemented  by  a  ATp -point  IFFT.  With  the  time-domain  channel 
estimate  Hls,  we  obtain  the  frequency-domain  estimates  using 
(17). 
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TABLE  I 

Input  Data  Structure  and  the  Corresponding  Bit  Rates 


K 

#  of  active 

subcarriers 

(Ka) 

#  of  null 

subcarriers 

(Kn) 

#  of  blocks 

in  a  packet 

(Nt) 

bit  rates 

without  coding 

2{Ka-Kl4) 

(T+Ta) 

bit  rates 

after  rate  2/3 

channel  coding 

512 

484 

28 

64 

10.52  kb/s 

7.0  kb/s 

1024 

968 

56 

32 

12.90  kb/s 

8.6  kb/s 

2048 

1936 

112 

16 

14.55  kb/s 

9.7  kb/s 

Stop  Packet  Stop  Packet  Stop  Packet  Stop 

,  K=512  ,  .  K=I024  ,  ,  K=2048  , 

N - H_ K  )l ^  ^ 


1^ -  3  packets  per  data  burst  - 

Fig.  3.  Each  data  burst  consists  of  three  packets,  with  K  —  512,  /C  =  1024,  and  K  —  2048,  respectively. 


4)  Multichannel  Combining:  Multichannel  reception  greatly 
improves  the  system  performance  through  diversity;  see,  e.g., 
[24]  on  multichannel  combining  for  single-carrier  transmissions 
over  UWA  channels.  In  an  OFDM  system,  multichannel  com¬ 
bining  can  be  easily  performed  on  each  subcarrier.  Suppose  that 
we  have  Nr  receive  elements,  and  let  and  de¬ 

note  the  output,  the  channel  frequency  response,  and  the  additive 
noise  observed  at  the  mth  subcarrier  of  the  rth  element.  We  thus 
have 


■  H^{m)  ■ 

'  vi,  ■ 

= 

d[m\  -h 

Nr 
_  '^m 

r'Nr 

L  m  J 

:—7.m  :=hm  Vm 

Assuming  that  has  independent  and  identically  distributed 
entries,  the  optimal  maximum-ratio  combining  (MRC)  yields 

d[m\  =  (h^hCj  h^Zm-  (21) 

Doppler  scaling  factor  CFO  and  channel  estimation  are  per¬ 
formed  independently  on  each  receiving  element  according  to 
the  procedure  described  in  Sections  III-B1-B3.  An  estimate  of 
the  channel  vector  is  then  formed,  and  used  to  obtain  the 
data  symbol  estimates  in  (21). 

IV.  Performance  Results  for  the  Experiment 
IN  Buzzards  Bay 

The  bandwidth  of  the  OFDM  signal  is  B  =  12  kHz,  and  the 
carrier  frequency  is  =  27  kHz.  The  transmitted  signal  thus 
occupies  the  frequency  band  between  21  and  33  kHz.  We  use 
ZP-OFDM  with  a  guard  interval  of  Tg  =  25  ms  per  OFDM 
block.  The  respective  numbers  of  subcarriers  used  in  the  exper¬ 
iment  are  K  =  512, 1024,  and  2048.  The  subcarrier  spacing 
is  A/  =  23.44,  11.72,  and  5.86  Hz,  and  the  OFDM  block  du¬ 
ration  is  T  =  1/A/  =  42.67,  85.33,  and  170.67  ms.  We  use 
rate  2/3  convolutional  coding,  obtained  by  puncturing  a  rate  1/2 
code  with  the  generator  polynomial  (23,35).  Coding  is  applied 
within  the  data  stream  for  each  OFDM  block.  QPSK  modula¬ 
tion  is  used.  For  K  =  512, 1024,  and  2048,  each  packet  con¬ 
tains  Nb  =  64, 32,  and  16  OFDM  blocks,  respectively.  The 


total  number  of  information  bits  per  packet  is  30  976.  The  signal 
parameters  and  the  corresponding  data  rates  are  summarized  in 
Table  I,  where  the  overhead  of  null  subcarriers  and  Kp  =  K/4 
pilot  subcarriers  is  accounted  for. 

Fig.  3  depicts  one  data  burst  that  consists  of  three  packets  with 
K  =  512,  K  =  1024,  and  K  =  2048,  respectively.  During  the 
experiments,  the  same  data  burst  was  transmitted  multiple  times 
while  the  transmitter  was  on  the  move. 

The  Woods  Hole  Oceanographic  Institution  (WHOI,  Woods 
Hole,  MA)  acoustic  communication  group  conducted  the  exper¬ 
iment  on  December  15,  2006,  in  Buzzards  Bay,  MA.  The  trans¬ 
mitter  was  located  at  a  depth  of  about  2.5  m  and  the  receiver 
consisted  of  a  four-element  vertical  array  of  length  0.5  m  sub¬ 
merged  at  a  depth  of  about  6  m.  The  transmitter  was  mounted 
on  the  arm  of  the  vessel  Mytilus,  and  the  receiver  array  was 
mounted  on  the  arm  of  the  vessel  Tioga.  OFDM  signals  were 
transmitted  while  Mytilus  was  moving  towards  Tioga,  starting 
at  600  m  away,  passing  by  Tioga,  and  ending  at  about  100  m 
away.  The  experiment  configuration  is  shown  in  Fig.  4. 

The  received  signal  was  directly  analog-to-digital  (A/D)  con¬ 
verted.  The  signal  received  on  one  element  is  shown  in  Fig.  5, 
which  contains  seven  data  bursts  or  21  packets.  The  following 
observations  can  be  made  from  Fig.  5. 

1)  The  received  power  is  increasing  before  packet  19,  and 
decreasing  thereafter.  This  observation  is  consistent  with 
the  fact  that  Mytilus  passed  Tioga  around  that  time. 

2)  A  sudden  increase  in  noise  shows  up  around  packet  19. 
This  noise  comes  from  the  Mytilus  when  it  was  very  close 
to  Tioga. 

3)  The  second  packet  was  severely  distorted.  The  reason  is 
unclear. 

Simple  data  processing  reveals  the  following. 

4)  The  signals  before  packet  19  were  compressed,  which 
agrees  with  the  fact  that  the  transmitter  was  moving  to¬ 
wards  the  receiver.  The  signals  after  that  were  dilated, 
confirming  the  fact  that  the  transmitter  was  moving  away 
from  the  receiver. 

Next,  we  present  numerical  results  based  on  the  sequence  of 
the  receiver  processing  shown  in  Fig.  1.  We  present  a  selected 
set  of  results  and  comparisons. 
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< -  600  m— 110  m - ► 


Fig.  4.  Configuration  of  the  experiment  in  Buzzards  Bay,  MA. 


Fig.  5.  Received  signal  (amplitude)  for  the  Buzzards  Bay,  MA,  experiment. 

A.  Doppler  Scaling  Factor  Estimation 

For  each  of  the  21  packets  transmitted,  the  algorithm  of 
Section  III-B 1  was  used  to  estimate  the  Doppler  scaling  factor. 
Based  on  each  Doppler  scaling  factor  a,  the  relative  speed 
between  the  transmitter  and  the  receiver  was  estimated  as 
V  =  a  X  c,  using  a  nominal  sound  speed  of  c  =  1500  m/s. 
The  relative  speed  and  the  resulting  Doppler  shift  at  the  carrier 
frequency  afc  are  shown  in  Fig.  6,  which  summarizes  the 
results  for  element  1 . 

We  see  from  Fig.  6  that  the  Doppler  shifts  are  much  larger 
than  the  OFDM  subcarrier  spacing.  For  example,  if  i)  =  8.30  kn 
(packet  15),  which  indicates  that  Mytilus  was  moving  toward 
Tioga  at  such  a  speed,  the  Doppler  shift  is  76.98  Hz  at  fc  = 
27  kHz,  while  the  subcarrier  spacing  is  only  A/  =  23.44, 1 1 .72, 
and  5.86  Hz  for  K  =  512, 1024,  and  2048,  respectively.  Hence, 
rescaling  the  waveform  (even  coarsely)  is  necessary  to  mitigate 
the  Doppler  effect  nonuniformly  in  the  frequency  domain. 

B.  High-Resolution  Residual  Doppler  Estimation 

The  high-resolution  CFO  estimation  was  performed  on  a 
block-by-block  basis,  as  detailed  in  Section  III-B 2.  Fig.  7 
shows  the  CFO  estimates  for  packets  5  and  17  for  iC  =  1024. 
We  observe  that  the  CFO  changes  from  block  to  block  roughly 
continuously  but  cannot  be  regarded  as  constant.  The  CFO 
estimate  is  on  the  order  of  half  of  the  subcarrier  spacing. 


Fig.  6.  Coarse  estimation  of  the  relative  speed  and  the  Doppler  shift  at  fc  — 
27  kHz  for  element  I . 


Fig .  7 .  Estimated  residual  Doppler  (CFO)  for  packet  5  (with  an  estimated  speed 
of  4.25  kn)  and  packet  17  (with  an  estimated  speed  of  8.26  kn).  The  CFO  fiuc- 
tuates  rapidly  from  one  block  to  another. 


Without  the  CFO  fine  tuning,  the  receiver  performance  would 
deteriorate  considerably. 

We  have  also  examined  joint  Doppler  scaling  factor  and  CFO 
fine  tuning  on  each  OFDM  block  based  on  null  subcarriers. 
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Fig.  8.  Channel  estimates  for  two  example  cases.  One  is  for  a  case  with  an  es¬ 
timated  speed  of  4.25  kn  (packet  5),  and  the  other  is  for  a  case  with  an  estimated 
speed  of  8.26  kn  (packet  17).  The  channel  delay  spread  is  about  4.5  ms.  There 
is  a  strong  direct  path  between  the  transmitter  and  the  receiver.  The  channel  en¬ 
ergy  in  the  8.26-kn  case  is  higher  than  that  in  the  4.25-kn  case,  as  the  transmitter 
is  closer.  The  second  peak  is  conjectured  to  be  from  the  bottom  bounce. 

which  requires  a  2-D  search  for  the  scale  h  and  the  CFO  e.  The 
performance  improvement  is  marginal  in  this  experiment,  so  we 
skip  the  results  on  the  joint  approach. 

C  Channel  Estimation 

Channel  estimation  is  based  on  equispaced  pilots,  as  detailed 
in  Section  III-B3.  Here,  we  use  Kp  =  K/A  pilot  subcarriers. 
Fig.  8  depicts  the  estimated  channel  impulse  responses  for 
two  cases.  In  one  case,  Mytilus  is  moving  toward  Tioga  at  an 
estimated  speed  of  4.25  kn  (packet  5),  and  in  the  other  case, 
it  is  moving  at  an  estimated  speed  of  8.26  kn  (packet  17). 
The  channel  duration  is  about  4.5  ms.  There  is  a  strong  direct 
path  between  the  transmitter  and  the  receiver.  The  energy  in 
the  8.26-kn  case  is  higher  than  that  in  the  4.25-kn  case.  This 
observation  matches  the  power  profile  shown  in  Fig.  5. 

A  second  path  is  also  observed  in  Fig.  8.  We  conjecture  that 
this  path  is  from  the  bottom  bounce.  This  conjecture  is  supported 
by  a  rough  computation  based  on  the  channel  geometry. 

Case  1)  Suppose  that  the  distance  is  400  m  and  the 
depth  is  12  m.  Then,  the  delay  between 
the  bottom  bounce  and  the  direct  path  is 
(2  X  ^2002  +  122  _  4oo)/1500  =  0.48  ms. 

Case  2)  Suppose  that  the  transmitter  is  now  150  m  from  the 
receiver  and  the  depth  is  12  m.  Then,  the  delay 
between  the  bottom  bounce  and  the  direct  path  is 
(2  X  ^752  +  122  _  150)/1500  =  1.3  ms. 

These  numbers  roughly  correspond  to  the  interarrival  times 
marked  in  Fig.  8.  The  arrival  corresponding  to  the  second  peak 
can  thus  be  assumed  to  be  from  a  bottom  bounce. 

D.  BER  Performance 

Now,  we  report  the  bit  error  rate  (BER)  performance  without 
and  with  coding.  The  Viterbi  algorithm  was  used  for  channel 
decoding. 


Fig.  9.  BERs  averaged  over  each  packet,  element  I.  Packets  10  and  19  (with 
K  —  512)  have  decoding  errors.  Packets  2,  14,  and  20  (with  K  —  1024)  have 
decoding  errors.  Packet  9  (with  K  —  2048)  has  decoding  errors. 


Fig.  10.  BERs  averaged  over  each  OFDM  block,  packet  19,  K  —  512, 
element  1. 


First,  we  plot  the  BER  averaged  over  each  packet  in  Fig.  9, 
for  one  receiver  (element  1).  In  total,  6  out  of  21  packets  have 
errors  after  channel  decoding.  We  now  look  into  the  BERs  for 
each  OFDM  block  inside  the  packets  with  decoding  errors.  The 
results  are  as  follows. 

•  Packet  2  has  22  out  of  32  blocks  in  error  after  decoding. 
This  received  packet  was  badly  distorted,  as  can  be  seen  in 
Fig.  5. 

•  Packet  9  has  4  out  of  16  blocks  {K  =  2048)  in  error  after 
decoding. 

•  Packet  10  has  2  out  of  64  blocks  (AT  =  512)  in  error  after 
decoding. 

•  Packets  14  and  20  have  5  out  of  32  block  (AT  =  1024)  in 
error  each,  after  decoding.  Except  packet  20  having  four 
consecutive  blocks  in  error  at  the  end,  the  error  blocks  for 
other  packets  are  sporadic. 
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Fig.  11.  Configuration  for  the  experiment  in  Woods  Hole  Harbor,  MA. 


•  It  is  interesting  to  look  at  packet  19,  which  has  17  out  of  64 
blocks  in  error  after  decoding.  The  BERs  on  the  block  level 
are  shown  in  Fig.  10.  The  major  portion  of  the  error  blocks 
occurred  when  the  transmitter  was  passing  by  the  receiver. 
As  we  observe  from  Fig.  5,  the  Doppler  frequencies  were 
changing  from  positive  to  negative  values  around  packet 
19,  and  the  noise  level  increased  considerably  during  the 
passing. 

We  emphasize  that  with  block-by-block  processing,  decoding 
errors  in  previous  blocks  have  no  impact  on  future  blocks,  as 
confirmed  by  Fig.  10. 

We  now  report  on  the  BER  performance  with  two  receivers 
(using  elements  1  and  2).  In  total,  there  are  four  packets  in  error 
as  follows:  packet  2  has  17  out  of  32  blocks  in  error,  packet 
9  has  1  out  of  16  blocks  in  error,  packet  19  has  14  out  of  64 
blocks  in  error,  and  packet  20  has  4  out  of  32  blocks  in  error. 
The  sporadic  block  errors  with  single-receiver  processing  are 
mostly  corrected  with  two-receiver  processing.  The  BER  plots 
are  omitted  due  to  space  limitations. 

For  a  real  system,  the  block  errors  could  be  corrected  via 
autorepeat  request  (ARQ)  procedures,  or  via  coding  strategies 
such  as  rateless  coding  [25,  Ch.  50]  that  can  effectively  handle 
lost  blocks. 

V.  Performance  Results  eor  the  Experiment 
IN  Woods  Hole  Harbor 

This  experiment  was  conducted  on  December  1,  2006.  The 
same  signal  set  as  described  in  Section  IV  was  used.  The  signal 
was  transmitted  from  a  depth  of  about  2.5  m  and  received  by 
a  four-element  vertical  array  with  interelement  spacing  0.5  m, 
submerged  at  a  depth  of  about  6  m.  The  transmitter  was  mounted 
on  the  arm  of  the  Mytilus,  and  the  receiver  array  was  attached 
to  a  buoy  close  to  the  dock.  OFDM  signals  were  transmitted 
while  Mytilus  was  moving  away  from  the  dock  starting  from 
a  distance  of  50  m  and  ending  at  about  800  m.  Then,  Mytilus 
moved  towards  the  dock.  The  configuration  is  shown  in  Fig.  1 1 . 

The  channel  condition  was  very  difficult  with  strong  multi- 
path  after  the  guard  interval  of  25  ms.  The  last  strong  path  is  ev¬ 
ident  at  about  80  ms,  as  shown  in  Fig.  12.  This  long  delay  spread 
is  likely  due  to  the  reflections  off  the  pilings  near  the  dock. 

With  the  channel  delay  spread  longer  than  the  guard  in¬ 
terval,  interblock  interference  (IB I)  emerges.  We  do  not  try  the 
channel  shortening  approach  to  reduce  the  IBI  before  OFDM 


LFM  correlation  results  for  the  Woods  Hole  Harbor  experiment 


ms 


Fig.  12.  Channel  response  estimates  obtained  by  the  linear  frequency-modu¬ 
lated  (LFM)  preamble  matching.  The  channel  in  the  Woods  Hole  Harbor,  MA, 
experiment  has  strong  returns  even  after  the  guard  interval  of  25  ms.  As  a  result, 
IBI  exists.  Unlike  this  situation,  the  channel  in  the  Buzzards  Bay,  MA,  experi¬ 
ment  has  delay  spread  much  shorter  than  the  guard  interval. 


demodulation  (e.g.,  using  methods  from  [26]-[28]).  Instead, 
we  treat  all  multipath  returns  after  the  guard  interval  as  additive 
noise;  hence,  the  system  is  operating  at  low  signal-to-noise  ratio 
(SNR).  Nevertheless,  with  channel  coding  and  multichannel 
reception,  reasonable  performance  is  still  achieved,  which 
speaks  for  the  robustness  of  the  receiver. 

To  illustrate  the  performance,  we  present  results  of  two 
data  bursts.  One  data  burst  was  transmitted  when  Mytilus  was 
moving  away  from  the  dock  at  a  low  speed  of  about  3  kn.  The 
other  data  burst  was  transmitted  when  Mytilus  was  moving 
towards  the  dock  at  a  high  speed  of  about  10  kn. 

A.  Doppler  Scaling  Factor  Estimation 

Table  II  shows  the  estimated  speeds,  which  reflect  the  exper¬ 
imental  settings.  The  Doppler  shifts  at  fc  =  27  kHz  are  very 
large  for  both  cases.  In  the  low-speed  case,  the  Doppler  shift  is 
on  the  order  of  the  OFDM  subcarrier  spacing  (23.44  Hz  when 
K  =  512).  In  the  high-speed  case,  the  Doppler  shift  is  much 
larger  than  the  subcarrier  spacing. 
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TABLE  II 

Coarse  Estimation  of  Doppler  Shift  and  Relative  Speed  for  Element  I 


The  low- speed  ease 

The  high-speed  ease 

Paeket 

Doppler  shift  due  to 
sealing  at  /c  (Hz) 

Relative  speed 
(kn) 

Paeket 

Doppler  shift  due  to 
sealing  at  /c  (Hz) 

Relative  speed 
(kn) 

1  (K=512) 

-23.84 

-2.56 

1  (K=512) 

91.49 

9.86 

2  (K=1024) 

-21.30 

-2.29 

2  (K=1024) 

87.88 

9.47 

3  (K=2048) 

-24.06 

-2.60 

3  (K=2048) 

96.03 

10.36 

Fig.  13.  Estimated  residual  Doppler  shift  of  packet  1.  K  =  512  and  each 
packet  has  64  OFDM  blocks. 


Fig.  14.  Estimated  residual  Doppler  of  packet  2.  K  =  1024  and  each  packet 
has  32  OFDM  blocks. 

B.  High-Resolution  Residual  Doppler  Estimation 

Figs.  13,  14,  and  15  show  the  CFO  estimates  for  packets  1, 
2,  and  3  of  element  1,  respectively.  The  following  observations 
are  made. 

1)  The  CFO  changes  from  block  to  block  smoothly,  but 
cannot  be  regarded  as  constant. 


Fig.  15.  Estimated  residual  Doppler  of  packet  3.  /C  =  2048  and  each  packet 
has  16  OFDM  blocks. 


2)  The  residual  CFO  effect  cannot  be  neglected. 

3)  The  CFO  estimates  are  on  the  order  of  half  of  the  subcarrier 
spacings  for  the  low- speed  case. 

4)  In  the  low-speed  case,  the  CFO  changes  periodically  over 
time.  The  period  is  the  same  for  all  three  settings.  In  the 
high-speed  case,  this  phenomenon  is  not  present.  A  pos¬ 
sible  explanation  for  this  effect  is  that  Mytilus  rises  and 
falls  due  to  waves,  which  is  more  pronounced  at  low  speed 
than  at  high  speed. 

5)  Note  that  fewer  null  subcarriers  are  available  in  the  K  = 
512  case  than  in  the  K  =  1024  and  K  =  2048  cases,  and 
hence  the  CFO  estimation  is  more  affected  by  the  noise  re¬ 
alizations.  When  K  increases,  more  null  subcarriers  lead  to 
better  noise  averaging,  and  the  corresponding  curves  look 
smoother.  This  trend  is  clearly  shown  in  Figs.  13-15. 

C.  Channel  Estimation 

Figs.  16  and  17  depict  the  channel  estimates  for  the  3-  and 
the  10-kn  cases,  respectively.  We  observe  several  stable  paths 
whose  delays  do  not  depend  on  the  location  and  the  speed  of 
the  transmitter.  For  example,  there  is  one  stable  path  around 
3  ms.  This  path  could  be  best  interpreted  as  the  first  reflected 
path  from  the  dock.  The  receiver  is  about  2  m  from  the  dock. 
Hence,  the  dock-reflected  path  will  be  delayed  by  2  x  2/1500  = 
2.6  ms  relative  to  the  direct  path.  This  is  a  constant  delay,  which 
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Fig.  16.  Estimated  channel  impulse  responses  (magnitude)  for  packets  1-3;  p;  jg  bERs  for  each  OFDM  block,  the  low-speed  case,  A'  =  2048. 
element  I  in  the  low-speed  case. 


Fig.  17.  Estimated  channel  impulse  responses  (magnitude)  for  packets  1-3; 
element  I  in  the  high-speed  case. 


Fig.  19.  BERs  for  each  OFDM  block,  the  high-speed  case,  K  —  2048. 


does  not  depend  on  the  distance  between  the  transmitter  and  the 
receiver. 

D.  BER  Performance 

Because  the  channel  condition  was  particularly  severe  in  this 
test,  both  coding  (rate  2/3)  and  multichannel  combining  were 
necessary  to  improve  the  BER  performance.  The  following  per¬ 
formance  results  are  obtained  with  three  receiving  elements. 

For  packet  3  with  K  =  2048,  Figs.  18  and  19  compare  the 
uncoded  performance  and  the  coded  performance  on  the  OFDM 
block  level,  with  single-channel  or  multichannel  reception,  in 
different  settings.  With  MRC,  the  uncoded  BERs  averaged  over 
the  packet  are  2  x  10“^  and  1.7  x  10“^  for  the  low-speed  and 
high-speed  cases,  respectively.  After  rate  2/3  coding,  the  BERs 
averaged  over  the  packet  are  1.6  x  10“^  and  5.8  x  10“^  for  the 
low-speed  and  high-speed  cases,  respectively.  We  observe  the 
following  from  Figs.  18  and  19. 

1)  The  uncoded  BER  is  large,  on  the  order  of  10“^  for  single¬ 
element  reception  and  10“^  for  multichannel  reception. 


2)  For  single-element  reception  with  large  uncoded  BER, 
coding  does  not  help.  However,  for  multichannel  recep¬ 
tion,  the  BER  performance  is  much  improved  when  coding 
is  used. 

With  K  =  1024,  the  BERs  averaged  over  the  packet  (packet 
2)  after  MRC  and  coding  are  1.1  x  10“^  and  6.5  x  10“^  for 
the  low-speed  and  high-speed  cases,  respectively.  With  K  = 
512,  the  BER  averaged  over  the  packet  (packet  1)  after  MRC 
and  coding  is  3  x  10“^  for  the  low-speed  case,  while  the  re¬ 
ceiver  does  not  work  well  for  the  high-speed  case.  These  re¬ 
sults  show  that  the  setting  with  larger  K  has  better  performance 
in  this  experiment.  When  K  increases,  the  effect  of  channel 
variation  within  one  OFDM  block  becomes  more  severe,  while 
on  the  other  hand,  the  receiver  has  more  null  subcarriers  and 
pilot  subcarriers  for  better  CFO  and  channel  estimation  against 
noise  (c.f..  Table  I).  Note  that  the  sampling  rate  is  fixed  for  all 
three  cases,  and  hence,  the  discrete-time  channel  has  approxi¬ 
mately  the  same  number  of  taps.  The  noise  effect  outweighs  the 
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channel-variation  effect  in  this  data  set,  since  the  receiver  oper¬ 
ates  at  a  noise-limited  region,  due  to  the  large  noise  contributed 
by  the  arrivals  after  the  guard  interval. 

Although  the  results  for  the  Woods  Hole  Harbor,  MA,  exper¬ 
iment  are  worse  than  those  for  the  Buzzards  Bay,  MA,  experi¬ 
ment,  they  demonstrate  the  robustness  of  the  proposed  receiver 
in  the  presence  of  a  difficult  channel  with  a  delay  spread  much 
longer  than  the  OFDM  guard  interval.  Note  that  a  16-state  rate 
2/3  code  is  used  here.  A  much  stronger  channel  code  (e.g.,  the 
nonbinary  low-density  parity-check  (LDPC)  code  used  in  129]) 
would  considerably  improve  the  BER  performance. 

VI.  Conclusion 

In  this  paper,  we  investigated  the  application  of  OFDM  in 
wideband  UWA  channels  with  nonuniform  Doppler  shifts.  To 
compensate  for  the  nonuniform  Doppler  distortion,  a  two-step 
approach  was  used:  resampling  followed  by  high-resolution  uni¬ 
form  compensation  of  the  residual  Doppler.  Null  subcarriers  fa¬ 
cilitate  Doppler  compensation,  and  pilot  subcarriers  are  used 
for  channel  estimation.  The  receiver  is  based  on  block-by -block 
processing,  and  hence,  it  is  suitable  for  fast- varying  channels. 

The  method  proposed  was  tested  in  two  shallow-water  ex¬ 
periments.  Over  a  bandwidth  of  12  kHz,  the  data  rates  were 
7.0,  8.6,  9.7  kb/s  with  QPSK  modulation  and  rate  2/3  convolu¬ 
tional  coding,  when  the  numbers  of  subcarriers  were  512,  1024, 
and  2048,  respectively.  Good  performance  was  achieved  even 
when  the  transmitter  and  the  receiver  were  moving  at  a  rela¬ 
tive  speed  of  up  to  10  kn,  where  the  Doppler  shifts  were  greater 
than  the  OFDM  subcarrier  spacing.  Experimental  results  sug¬ 
gest  that  OFDM  is  a  viable  candidate  for  high-rate  transmission 
over  UWA  channels. 

Future  research  will  address  several  topics,  including  short¬ 
ening  methods  for  channels  whose  delay  spread  is  longer  than 
the  guard  interval,  extension  of  resampling  to  generalized 
time-varying  filtering  for  channels  with  different  Doppler 
scaling  factors  on  different  paths,  and  multiple-input-mul- 
tiple-output  (MIMO)  techniques  129]-131]. 
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